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ABSTRACT. Two major deviations from causality in the existing formulations of 
quantum mechanics, related respectively to quantum chaos and indeterminate wave 
reduction, are eliminated within the new, universal concept of dynamic complexity. The 
a analysis involves a new paradigm for description of a system with interaction, the 
^ principle of dynamic multivaluedness (redundance), and the ensuing concept of the 
^ fundamental dynamic uncertainty. It is shown that both the wave reduction and truly 
unpredictable (chaotic) behaviour in quantum systems can be completely and causally 
^ understood as a higher sublevel of the same dynamic complexity that provides the 
^ causally complete picture of the unified wave-particle duality and relativity at its lowest 
2 level (quant-ph/9902015,16). The presentation is divided into five parts. The first three 
^ parts deal with intrinsic randomness in Hamiltonian (isolated) quantum systems as the 
^ basic case of dynamical chaos. In the last two parts a causal solution to the problem of 
quantum indeterminacy and wave reduction is proposed. Part I introduces the method of 
0-{ effective dynamical functions as a generalisation of the optical potential formalism. The 
g method provides a legal transformation of the Schrodinger equation revealing the hidden 
^ multivaluedness of interaction process, i. e. its self-consistent, dynamical splitting into 
^ many equally real, but mutually incompatible branches, called 'realisations'. Each 
realisation incorporates the usual "complete" set of eigenfunctions and eigenvalues for 
the entire problem. The method is presented in detail for the Hamiltonian system with 
periodic (not small) perturbation, in both its time-independent and time-dependent 
versions. 



NOTE ON NUMERATION OF ITEMS. We use the unified system of consecutive numbers for formulas, 
sections, and figures (but not for literature references) throughout the full work, Parts I-V. If a reference to 
an item is made outside its "home" Part of the work, the Roman number of this home part is added to the 
consecutive number: 'section 9.2' and 'section 9.2.IV' refer to the same, uniquely defined section, but in 
the second case we know in addition that it can be found in Part IV of the work. 



* Address for correspondence: Post Box 1 15, Kiev - 30, Ukraine 252030. 
E-mail address: kiril@metfiz.freenet.kiev.ua 



1. Introduction 



The recent emergence and development of the dynamical chaos concept have 
engendered profound changes in our understanding of both fundamental and 
practical aspects of the dynamical system behaviour in many different fields of 
physics (e. g. [1]). In particular, the behaviour of simple non-dissipative 
mechanical systems with few degrees of freedom presents an elementary case 
well suited for the study of the fundamental origins of dynamical randomness, 
with possible further extension to more complex situations. The description of 
chaos in such elementary dynamical systems within the formalism of classical 
mechanics has seemed to be rather successful and self-consistent [1-4]. At the 
same time its proposed quantum-mechanical versions, despite a large amount of 
the efforts made, have failed in creating a similar prosperous situation, even 
though a number of important particular results has been obtained [4-7]. The 
problem of the very existence of the truly unpredictable behaviour of 
deterministic quantum systems remains unsettled [8]. The fundamental 
difficulty stems from the unavoidable wave involvement in quantum postulates: 
waves do not easily show global instabilities necessary for the development of 
chaotic regimes or, in other terms, waves lead to discreteness (one may 
physically realise minimum a half-wave and not, say, one fifth of it), and the 
discreteness is incompatible with the existing notion of instability appealing to 
infinitesimal values and apparently indispensable for the known definitions of 
chaos. Solutions to this basic problem, proposed or implied, vary over a wide but 
already shrinking range. 

The most popular point of view involves the reduction of quantum chaos to 
some kind of very intricate but basically regular behaviour related eventually to 
the peculiarities of involved mathematical objects like zeta-functions [4,9] or 
random matrices [10]. Whatever the variations of such approach, it suffers from 
the evident conflict with the rather attractive idea, strongly supported by the 
existence of different types of classical chaos, that the world's complexity 
should be greater than zero: the latter would be impossible without the 
irreducible 'true' unpredictability and randomness in quantum mechanics as the 
most general description of the world. If it is not the case, i. e. if, for example, 
algorithmic complexity of quantum dynamics turns invariably to zero, then one 
possibility is that the existing quantum mechanics is no more general enough to 
provide a non-contradictory description for the world of chance [11]. 

The practical absence, however, of such properly chaotic scheme of 
quantum mechanics leaves the way open for another possibility: the conventional 
quantum mechanics itself is valid, for classically chaotic quantum systems in 
the same sense that for the classically regular ones, while the chaoticity of its 
classical limit "appears" due to the particularities of the semiclassical 
transition. The emerging "quantum chaology" [12] finally also proposes rather a 
possibility than a solution: it is suggested why quantum and classical 
descriptions need not provide similar concepts of chaos, but it remains unclear 
why and how exactly the quantum-mechanical reality with complexity zero 
(absolute dynamical predictability and reversibility) passes to its truly random 
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limiting case witli non-zero complexity at the other end of the semiclassical 
transition. The latter should therefore 'produce' complexity 'from nothing' by a 
still unknown mechanism. 

This gives us also a hint to the third and logically the last remaining 
possibility: the true solution is hidden neither within the semiclassical 
transition, nor in a reinterpretation of quantum mechanics, but rather in our 
understanding of chaos and complexity in general, even in classical mechanics. It 
may imply that in reality the 'pure' randomness and the corresponding 
complexity do not exist at all, and the apparent world of chance is just an 
illusory product of the ever wavering spirit... This choice, rather exotic as it may 
seem, can be provided, nonetheless, with a precise rational basis: it is stated 
simply that as the quantum-mechanical "diffusion localisation time" [13], after 
which a system starts deviating from apparently chaotic behaviour, becomes 
practically infinitely large for classical objects, one can never have a chance to 
achieve this regime of chaos suppression in the classical world, contrary to 
small quantum objects (this opinion was communicated to the author by Dr. S. 
Weigert from Basel University). The world is thus not complex, within this 
possibility, neither quantum-mechanically, nor classically, even though it is 
permitted to be very intricate. It represents, in fact, a huge generator of random 
numbers, basically deterministic but very perfect in its simulation of 
randomness (in other words, the world's dynamics is a periodic one, but the 
period is extremely large). In this way the problem of the origin of quantum 
chaos is effectively reduced to the alternative between the true complexity and 
a simple intricacy of being. 

Apart from this logical choice among the three basic possibilities, one may 
mention the hypotheses about the quantum measurement involvement in chaos 
(see e. g. the special NATO Workshop Proceedings [14]) or the "fundamentally 
irreducible" representation by the density matrix [15]. It is not difficult to see 
that these approaches can also be reduced to one of the above possibilities 
touching the roots of the quantum paradigm. At their present state, however, 
they seem to be either too abstract, or too special to form a universal basis for 
the first-principle understanding of quantum chaos in real physical systems.*) 

In the absence of clearness in the fundamental aspects, there is the general 
tendency to study particular features, or "signatures", of chaos in the regular 
quantum dynamics of the classically chaotic systems [6,7] leaving aside the 
related logical puzzles. However, it becomes more and more evident that one 
cannot hope to be really successful even in such particular research without the 



*) We do not discuss here the appearing multiple attempts to relate deterministic randomness in 
quantum systems to a specific behaviour of the dissipative quantum systems inevitably involving an 
external source of irreversibility, whether it is stochasticity of the environment or the equivalent 
"coarse-graining" of a problem. We are looking for an intrinsic causal source of dynamic 
unpredictability which would remain in the zero-noise limit, even though the addition of noise may play 
an important practical role in complex system behaviour (see sections 2. 3. II and 5. Ill below). This 
approach involves eventually considerable extension of the existing interpretation of any type of chaotic 
behaviour (see sections 6. Ill and 10.V of this work), including the "divergent-trajectories" paradigm 
of classical deterministic chaos that suffers, in fact, from the same incompleteness as the "dissipative 
chaos" concept. 
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general self-consistent picture including the basic issues (not to mention the 
temptations of pure curiosity). This situation predetermines the importance of 
the search for a general self-consistent quantum (and eventually any other) 
chaos description providing the fundamental origin of randomness in 
deterministic systems. In this essay we present such an approach starting from 
the application of the unreduced version of the well-known optical potential 
method (see e. g. [16]) to the analysis of quantum chaos in Hamiltonian systems 
with periodic perturbation. This approach appeared originally as a part of 
quantum-mechanical description of charged particle scattering in crystals [17] 
revealing its chaotic behaviour (see, in particular, section 2.5 of the cited 
article). The results obtained are considerably developed and generalised in the 
present paper. 

Using a generic example of arbitrary periodically perturbed Hamiltonian 
system, we show that our method naturally leads to the new concept of the 
fundamental origin of chaos in dynamical systems (section 2). In particular, it 
permits one to overcome the 'pathological regularity' of quantum mechanics and 
to perform the ordinary semiclassical transition also for chaotic systems, in 
agreement with the correspondence principle (section 3. II). These results show 
that our approach is presented in the form ready for its application to practical 
study of Hamiltonian quantum chaos in various real physical systems including 
both basic aspects and the particular analysis of the measured quantities. 

An important topic concerning fractals involvement with chaos is 
considered in section 4. Ill, where we show, within the same approach, that 
fractals naturally appear as solutions of the modified Schrodinger equation for a 
chaotic system that can be obtained analytically and not only by computer 
simulations (the result refers, in fact, to any dynamic equation presented in the 
proper form). In fact, we demonstrate that any solution describing a system in a 
chaotic regime (and this is the absolute majority of all situations) has a fractal 
character, and the properties of this fractal can be determined by the outlined 
procedure. 

It is important that in order to obtain quantum dynamics with non-zero 
complexity and the conventional semiclassical transition, one does not need to 
reconsider the foundations of quantum mechanics as such but rather to use 
another, more general, form of the same formalism. This can help to moderate 
the painful choice described above (we discuss these issues in more detail in 
section 5. III). Moreover, in parts IV, V of this work we show that it is the main 
unsolved problems of the foundations of quantum mechanics, known as quantum 
indeterminacy and wave reduction, that can be given transparent causal 
solutions by application of the same method to the process of quantum 
measurement. 

Finally, the purpose of this work is to demonstrate the universal character 
of the results obtained permitting one to extend the same concept of complex 
behaviour to other Hamiltonian and non-Hamiltonian, classical, and eventually 
distributed nonlinear systems. The ensuing universal notions of the fundamental 
dynamic uncertainty, randomness, probability, complexity, (non)integrability, 
and general solution are introduced in section 6. III. 
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2. Formulation of the method 
2.1. Effective dynamical functions 



Consider a conservative dynamical system witli tine Hamiltonian 

H = h + V, 

wliere h is tlie free motion Hamiltonian, and V describes the elastic interaction. 
Within the total Hamiltonian H we separate the integrable part, H q, 
corresponding to the regular dynamics, and the perturbation, Hp, (generally, not 
small) inducing chaotic behaviour of the whole system which we want to 
describe: 

H=Ho + Hp, Ho^ho+Vo, Hp^hp+Vp . 

We introduce then a particular representation, for definitness chosen in the 
form of coordinate representation, H = H(r), and divide the vector of independent 
variables into two parts, r = {ra.r^}, so that, in accord with the integrability of 
Hq, Ho = Hoir^) and Hp = Hp{ra,rTi). If chaos is induced by the addition of extra 
degrees of freedom (dimensions), the variables r-^ may correspond to these 
degrees, while the motion limited to the degrees Is considered to be regular. 
In cases where such subdivision is not naturally given by the conditions of a 
problem, it can always be made using the well-known regularity of one- 
dimensional problems: in the simplest version a one-dimensional component of r 
is chosen as r^, Va = x, r = {x,rn}. The subdivision of the Hamiltonian is then 
performed with the help of the Fourier analysis or other suitable expansion. As 
we shall see later, the case of time-dependent perturbation can also be 
considered within the same formalism, and then rjt corresponds to the time 
variable, t. If the degrees of freedom remain unchanged and chaos is due to a 
symmetry-breaking perturbation, one can still use the above formal method of 
division. To obtain physically more meaningful result, one may imply, where 
possible, the quantum analogues of the action-angle variables (r^ I, r-^ ^ 0), 
or those determined by the symmetry (e. g. = r), or other suitable choice.*) 

While studying purely dynamical origins of stochasticity, it is natural to 
assume that Vp is a periodic function of r-^. It is also one of the most interesting 
particular cases (e. g. particle scattering in regular structures [17,18], atom 
excitation by electromagnetic radiation [19], the kicked rotor model [20], etc.), 
and it is well suited for the demonstration of our method. The dependences of Vo 
and Vp on may be periodic or not. 



*) We consider it to be always possible. Although our description encompasses, in principle, this latter 
type of chaos, the results below are specified rather for the former one, with some generalising notes 
where possible; the detailed investigation of the second type of chaos is left for next publications. 
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Consider now the Schrodinger equation for our system, 

[hoira) + ^p(r„) + Voir a) + Vplr^.r^)] "f^ir^.r^) = EW{ra,r^) , (1) 

where E is the total energy. The involvement of periodicity inspires the idea 

about the Fourier transformation over r^*) which is done in the usual way and 
leads to the system of equations for the component functions \j/g^{ra): 

[hoira) + Vo{ra)]Wo{ra) = £aWo{ra) - T^V-gni^o)WgT,ira) , (2a) 

^o(ra)V^gJra) + Y.^g^-Sn'i^o)Vgni^o) = Eog^WgJj^o) " Jr^) V^o(ra) , (2b) 

where g^, g^t' ^ are the dual "reciprocal lattice" vectors with respect to the 
"direct lattice" of vectors r^, 

£ag^ = E - h2(K;, + gn?l2m , £a ^ E - Yi^K^^Ilm , (3) 



W{r) = exp(iK;trn) [v/-o(ra) + V/-g^(ra)exp(ignrn)] 
V{r) = Vo(ra) + l^Vg;t(ra)exp(ignrn) , 



and the wave vector corresponds to the standard Bloch representation of the 
wave function for periodic potentials (see e. g. [16]). 

Now one can use this decomposition to study the influence of "chaos 
bringing" perturbation Vp(ro,r;t) providing the terms with g;t ^ 0. We start by 
applying the simple method of substitution and first express YgA^o) through 

y/o{ra) from eq. (2b) with the help of the Green function for its homogeneous part 
with respect to yfg^{rc), 

hoiroWgni^o) + I^Vgn-gni^oWgni^o) = £ag^¥gni^o) ■ (4a) 

The Green function is given by the well-known expression: 



I 



0* 

¥gnniro)¥^Jra ) 



^g3t(^0'''o ) - ^ J 

%7t« ■ ^(^gn 



n 



*) In the general case it will be an expansion in terms of the other complete system of functions 
appropriate to a problem. 
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where {v^g^„(ra)} and {Cg^jw} are the sets of eigenfunctions and eigenvalues, 

respectively, for the auxiliary system of equations (4). The solution of the 
system (2b) can be expressed as 

Vgn(^o) = - j^ra'Gg^(r(j,r(j')Vg„(r(j')V^o(ra') . 

So 

where the domain of integration coincides with the "unit cell" for V{r) 
periodic in or with the whole domain of definition on re? for a non-periodic 
potential. 

Now we substitute the obtained expression for Wgni^o) in the right-hand 
side of eq. (2a) and come to the conclusion that the problem is reduced to the 
modified Sclirddinger equation for v^o(ra): 

[/Zo(ra) + Veff(ra)]V^o(ra) = eaV^o(ra) , (5) 

where the ordinary potential Vo(ra) is replaced by the effective potential (EP) 
Veffiva), also known as coherent, or optical, potential [16]. It is obtained as a 
sum, 

Veff(ra) = Vo(ra) + i^{ra) , (6a) 
where i}{ra) is the nonlocal part of EP expressed by the integral operator: 

7^(ra)/(ra)- I Jra'y(ra,ra')/(ra') . (6b) 

So 

The integral kernel V{ra,ro') can be presented in the form 



n* 

V{ra,ra) = / J , , (6c) 

- eg;t« - ejig;t - 2cos«g^'y (£ - ea)e7tg;t 



§71 ■ 

where 



2 2 



and ttg^ is the angle between the vectors and (for one-dimensional r^, a^^ 

takes only two values, a^^ = 0,n). 

This expression for EP contains, in particular, the unknown solutions of the 
auxiliary system of equations (4). Their properties are investigated in Appendix 
(see also [17]), where it is shown that qualitatively they are similar to those of 
solutions of the Schrodinger equation with the unperturbed potential Vo(ra). More 
detailed quantitative characteristics of these solutions and their properties can 
be analysed using suitable approximations, but we need not use much of these 
details here because our conclusions, as we shall see, depend little on them, but 



7 



rather on the general form of representation (5)-(6) of a problem. In particular, 
the effective dependence of these unknown solutions on is rather weak, for 
the significant terms of the sum over in eq. (6c), which permits us not to 
specify this dependence (see also Appendix). 

To obtain the complete solution to a problem, one should find the solutions 
of the modified Schrodinger equation, eq. (5), and then substitute them into the 
expression for y/g^irc), after which the general solution can be written as 

•F(r) = \ Cn[yfon{ra) + 12 V^g;tn(ra)exp(ig;tr;t)]exp(iK;t„r;t) = 
n 

= / c„exp(iK„„rn)[i + ^exp(ignrn)(§g^„(r(j)] V^on(ra) , (7) 



where 



¥gnni^c) = ^g;in(ra)V^on(ra) = \ draXg^n{rc,rc)¥o{ro') 



y 



(8) 



the coefficients c„ are to be determined from the boundary or initial conditions, 
as well as certain components of K, £g^n and Kjt„ are specified by eqs. (3) with Sa 

= £an, and {y/on}, {Son} ^fo the complete sets of eigenfunctions and eigenvalues for 
the modified Schrodinger equation, eq. (5). And finally the main measured 
quantity, the probability density distribution (PDD), p(r) = |'F(r)p, can be 
expressed directly from eqs. (7), (8). The resulting general formula for it can be 
found elsewhere (see eq. (18) in ref. [17]). 

Before analysing the results obtained, it would be not out of place to note 
that another practically important case, that of time-dependent periodical 
perturbation. Hp = Vp{ra,t), is effectively described by the same system of 
equations (2), where one should make the substitutions 

t, k {k ^ \s an integer), Sag^ £ok = £o - ^(O-nk, 

the wave function being presented in the form 

W{ra,t) = exp(-iea?/h) [V^o(ra) + L Wk{ro) eKp{i COM] , 

k 

and (On being the main frequency of the perturbation: 
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k 

This conclusion can be verified starting from the substitution of the total wave 
function above into the time-dependent Schrodinger equation. Then it is easily 
seen that this problem can also be reduced to solution of the modified 
Schrodinger equation (5), where the kernel of the effective potential can be 
expressed in a slightly different form, 



I 



V.k{ra)Vk{ral¥kn{ra)wT{ral 
£c - £kn - hCOnk 



k,n 



and {Wkni^o)}, {^kn) s""® determined from the auxiliary system of equations: 



ho{ro)Vk{ro) + Y.^k.k'{roWk'{ra) = £okVk{ro) , Kk' ^ 0. (4b) 

k' 



The general solution is (cf. eq. (7)): 



'^'('■a.O = 7^ C„exp(-iean?/h)[l + 1^ exp(ift)nA;?)^A:n(ra)] V^on(ra) 



with the integral kernel xkn{Yo,ro') of the operators t^kni^o) given by 



I 



Xkn{Yo,Yo') = Vkird / J , 

n ' ^kn ' ' h ft) 



and {£an}, {V^on(ra)} determined from eq. (5). Starting from these formulas one 
obtains any desired measurable quantity, similar to the time-independent 
formalism. 

It is clear that the time-dependent case is a generalisation of the well- 
known group of the model "kicked" systems with similar behaviour described by 
the basic standard map [20] (called also the standard model, see e. g. [2,3]). This 
model corresponds to the 5-like periodic kicks (i. e. the components Vk{ra) = 
Vo(ra) do not depend on k) and also to some particular choices of the potential 
dependence on = x (typically Vo{x), Vk{x) oc -cos(.x;)). The proposed 
generalisation can serve thus as a more realistic representation for many 
particular physical systems studied with the standard model. By analogy one may 
designate our time-dependent case as the generalised kicked oscillator. Because 
of its similarity to the time-independent case we now continue the analysis of 
the latter mentioning the differences between the two where necessary. 
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2.2. Fundamental dynamic multivaluedness 



It is not surprising to see, from the above expressions, tliat tlie EP metliod 
thus formulated cannot provide directly the exact solutions and is nothing but 
another formulation of a problem. However, we can show now that for the 
chaotic systems it is this representation that is much more relevant than the 
ordinary one (in this case, the Schrodinger equation with the ordinary potential, 
V(r)). It permits one to obtain, in a natural and self-consistent manner, a basic 
source of randomness and complexity in such systems and then to study, in 
relation to these fundamental concepts, their particular chaotic properties in 
terms of observable quantities. 

We see from the above formulas that the experimentally measured PDD is 
determined by the dynamics in EP. The distinctive property of the latter, and the 
most important one as far as the dynamical chaos is concerned, is its self- 
consistent dependence on the energy eigenvalues to be determined. This 
dependence appears in the explicit form when one tries to find the energy 
eigenvalues {£an} from the modified Schrodinger equation, eq. (5), while using the 
expressions (6c, d) for the integral kernel of Veff(ra) which depends itself on £a- 
Note that this property is restricted solely to the full non-perturbative EP 
formalism studied here as opposed to its various perturbative versions 
eventually used in many applications [16] including the problem of quantum chaos 
[21,22] (see also section 5. III). In fact, we deal here with the intrinsic effective 
nonlinearity of a chaotic system which is not taken into account either by the 
ordinary, non-modified formalism, or by the perturbative approaches (further 
discussion of the effective nonlinearity can be found in sections 5. Ill, 6. Ill, and 
9. 2. IV). Now we are going to show that this peculiar property leads to the 
conclusion that, instead of eigenvalues and eigenfunctions for the Schrodinger 
equation with the "non-chaotic" potential Vo(ra). one obtains up to 

A^max= {NnNn + l)Na (9a) 

solutions for the equation with EP, where Nji and Nti' are the numbers of terms in 
the sums over gjt and n, respectively, in eq. (6c). Among these A^max solutions, A^o 
= (A^7t' + l)A^a solutions correspond to the normal set of eigenfunctions for the 
full-dimensional Schrodinger equation with the potential V(r) (in a standard 

situation there should be A^^ = N^', neglecting the eg^„ dependence on g^, see 
Appendix). The additional growth of the number of solutions by 

A^A - A^max - A^o = N^'iNn " > Na 

cannot be explained in terms of ordinary splitting effects. These extra solutions 
could be, in principle, spurious, unphysical ones (e. g. unstable). This can indeed 
happen for some ranges of parameters, and this is precisely one of the 
mechanisms of regularity in chaotic quantum systems, studied below. However, 
it is difficult to imagine, and we confirm these doubts below, that this may be 
the case for all parameter values. We shall see that most often at least some of 
these additional solutions are quite real and observable. 
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Before proving these statements, note that it will be convenient to count 
the solutions in terms of Nc-Nr = N/Na, where N is the total number of 
eigenvalues (eigenfunctions) in a general case; 

aC'''^ A^max/A^a = A^A' + 1; Nr = No/No = Nn +1; (9b) 

= N^^^ - N^r = Na/No = Nn{Nn - 1) > 1. One may introduce also Ur = A^r/N? and 

^max ^ ^max^^O ^ ^ + 1) ^ A^;, for Nn » 1 ■ (9c) 

The statements above, concerning the number of solutions, can be verified 
by at least three different ways giving all the same result, eqs. (9): one can 
count solutions directly by analysing the system of equations corresponding to 
eq. (5) in the momentum or other suitable representation; one can study a 
problem by a graphical method; at last, reasonable and easily treated 
approximations of eqs. (5)-(6) can be proposed. Consider now, in consecutive 
order, these arguments in detail. 

The momentum representation of eq. (5) is 

ga' 

where, for the case of V(r) periodic in r^, the notations h„ , , and V^^^ ,hold 

so so go, go 

respectively for the Fourier representations of the free-motion Hamiltonian, the 
periodic part of the wave function v^o(ra) and the EP kernel y(ra,ra'), 

^-gTt. go-go ^gTl.go +go "^gjlgo ^gTtgo ' 

go". go'" 

, (11) 

gn,n 

while the (straightforward) details of the definitions of the participating 
Fourier-components of the total potential and wave function are not important in 
the context. For the case of V(r) non periodic in one can use expansion in any 
appropriate complete set of orthogonal functions (for example, the set of 
eigenfunctions for Ho), which gives effectively the same result with g^, ga', ■■■ in 
eqs. (10), (11) enumerating these functions. Now if we neglect £a under the root 
in the denominator of eq. (11) for larger E, E » £a, then diagonalisation of system 
(10) is evidently reduced to solution of an algebraic equation in of degree A^max 
defined by eq. (9a). The latter statement is thus valid for the maximum 
attainable value of the master equation degree and the number of solutions of 
the modified Schrodinger equation, while the influence of the neglected 8a term 
is considered in detail below, in section 3. II. 
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In what follows we call this additional splitting of solutions and, 
correspondingly of the effective Hamiltonian, potential or other relevant 
dynamical function, the fundamental multivaluedness of dynamical functions 
(FMDF). It includes the ordinary dimensional splitting into the complete set of 

Nr solutions as one component, but also, possibly, other similar complete 

components. We shall call the i-th component the Hh realisation of a dynamical 
system (or of a problem) and designate it by 9tj; it effectively includes, besides 
the corresponding complete set of solutions, also the corresponding components 
of EP and PDD: 

- {{£cnV, {¥on{rc)V, V'eff(ra). P/W} {i = 1,2,.. .iVg^) , (12) 

where > 1 is the total number of realisations in their ensemble thus obtained. 
It is formally limited from above: 

N^i < N^^"" <Nr-N%l< nT^"" -N^, + 1 = Nn'iNn - 1) + 1 

(we shall see below that actually = A^^). 

Now we can present another, the most convincing and meaningful, evidence 
supporting the reality of the fundamental multivaluedness of solutions and based 
on the graphical analysis of eq. (5). We first rewrite eq. (5) for certain n-th 
eigenvalue: 

[/Zo(ra) + Vo(ra)]V^on(ra) + j^^ra'V„(ra,ra')V^o«(ra') = £a«V^o«(ra) ■ 



Multiplying it by Woni^o) and integrating over we arrive at the following 
formulation of a problem: 



where 
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Vnn{£on)= / . o , (14a) 



Vll = /^i?raV^g^„-(ra)Vg^(ra)V^o«(ra) , 



s 



and 



ea«= \ dra%Jra)[ho{ra) + Vo(ra)] V^o«(ra) 
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Function V„„(£a«) possesses a number of singularities determined by the zeros of 
the denominators of each term of the double sum in eq. (14a). The zero, and thus 
the corresponding singularity position, for a term with certain n' and gn is 

easily found to be at £an = £an'(gjt)' where 

(15a) 

The superscript "±" in the notations above serves to remind us about a feature 
important for the following analysis: for each g-^ there are two terms in the sum 
(14a) with the opposite directions of the vector corresponding to change of 
sign of costtg^ in eq. (15a). This can be expressed in a straightforward fashion 

for the simple and rather common case of one-dimensional perturbation (rn^z): 

For the case of time-dependent perturbation (the generalised kicked 
oscillator) one obtains, using the substitution procedure described above, the 
same eq. (13) with 




(14b) 



k,n ' 



where = Vg'^ \sn = k positions of singularities are determined as 

= ^kn- + "^COnk = Ekn ± hCOM . (15b) 



Now to apply our graphical analysis we plot, in Fig. 1 (p. 15), the left- and 
right-hand sides of eq. (13) vs £an taken as continuous independent variable. The 
solutions are found then as abscissae of the points of intersection of the curves 
corresponding to the two functions plotted. As is seen from the figure, 
representing eq. (13) for two characteristic cases of parameter values (see 
section 3. II), function V„„(£an) consists of many branches due to the sums over n' 
and gn of divergent terms. The points of divergence are designated by the 
corresponding vertical asymptotes, their positions on the horizontal axis being 
determined by eqs. (15). It is this plurality of branches which gives multiple 

solutions for e^n (we enumerate them by another index, j: ejj^, j > 0) instead of the 

single one, £an = ^aw for the Schrodinger equation with the 'regular' potential 
Vo{ro)- One part of this splitting can be explained by the ordinary multiplication 
of the number of eigenvalues and eigenfunctions due to the addition of the 
degrees of freedom corresponding to to those of Va- However, there is also 
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another part of the splitting which just represents the FMDF described above. 
Thus for the case of Fig. 1(a), where we have restricted ourselves only to two 
terms of each of the sums over n and ±\g%\ (i. e. Nti = 4, N-n' = 2), in the absence of 

this additional splitting one could not obtain more than A^^^A^^^' + l =3 

solutions, whereas actually one can count Nr = n'^^^ = 9 solutions for Son in full 

agreement with the general expression, eq. (9b).*) To deduce eq. (9) from our 
graphical analysis, we just note that, as follows from eqs. (14), (15), the number 
of asymptotes is equal to NjiN-n:', and then, as is clear from the figure, the number 

of points of intersection of the two curves, and thus of the solutions, is N^^^ = 
NnNn^ + 1- We shall see in section 3. Ill why and how this maximum number of 

solutions can diminish down to N^r giving rise to qualitative changes in the 

chaotic behaviour of a system. Figure 1(b) reproduces these results for another 
characteristic parameter values and for the value N^' = 3 (we continue this 
analysis in section 3. III). It can be especially clearly seen from configuration in 
this figure that the eigen-solutions are grouped so as to form N^r = A^^ 

realisations. Note that the constant adjusts itself, in a self-consistent 

manner, to its particular value for each realisation.*) 

It would be useful to introduce the notion of distance between two 
realisations. For our purposes it is quite sufficient to define it approximately as 
the average, or typical, magnitude of the difference between the values of the 
corresponding branches of EP. Then it follows from the above analysis that the 
separation between realisations (i. e. the distance between the neighbouring 
ones) is generally between the energy-level separation for the unperturbed 
potential, Ae^, at maximum (in the situation of global regularity. Fig. 1(b), see 

section 3. Ill), and the minimum of 2(£T^g^QE)^^^ = lizli^lElm Id^ (in the situation 

of global chaos. Fig. 1(a)). We see that the finite separation between the 
realisations, proportional to h, (i. e. their discreteness) is a quantum effect. As 

will be shown in the next section, it has the important physical consequences. It 
means also that in the semiclassical situation one has, at the scale of the 
characteristic potential values, many closely separated realisations forming a 
quasi-continuous distribution. Of course, these generic rules do not exclude the 
existence of particular realisation separations of greater (or smaller) 
magnitudes. 



*) Note that for convenience of illustration we consider cases with A^^ ^ N-n', whereas normally one 
should have A^^ = Nn'; however, this choice does not influence the obtained conclusions. 

*) As the value of may vary for different solutions e'cr^, the segments of the line San - ^on' 

intersecting the respective branches of the function y„„(ea„) in Fig. 1, can be slightly displaced 
vertically one relative to another. However, this will not produce any significant changes on the scale of 
our schematical representation, neither in the conclusions obtained, and we do not show these secondary 
details to avoid unnecessary complication. 
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Fig. 1. Graphical solution of the modified Schrodinger equation, eq. (13), for the Hamiltonian system 
with periodic perturbation in the parameter domains of global chaos (a), and global regularity (b). We 
plot the left- and the right-hand sides of eq. (13) vs eigenvalue to be determined, eo„. The illustration 
corresponds to the case E » ean (see text) and the following numbers of terms in the sums over and 
in the expression for V„„(eon), eq. (14a): = 4, Nt^' = 2 (a), and N,^ = 4, Nt^^ = 3 (b). The asymptotes 
of the function V„„(e(j„) are shown by the dashed lines with the respective values of g and n' marked at 
the lower end of each of them, where the non-zero integer g enumerates g^: = Ing/d^^ {g^^o = In/dT^). 



Two vertical dash-dotted lines in (a) correspond to the values eo„ = Sg^n'- 
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Finally, we refer to our paper [17] for the detailed description of a simple 
model for eq. (13) using another form of Vnni^an)- This model corresponds to the 
simplest choice A^^ = 2, A^^' = 1, A^^ = 1 and is reduced to a cubic equation that can 
be directly analysed including the dependence on parameters. Despite apparently 
rough approximations, the model is shown to be qualitatively reasonable for a 
large range of practical cases. We do not reproduce it here because effectively it 
represents a particular case of a more general analysis performed above. 
Nonetheless, it gives a transparent illustration of the appearance of extra 
solutions and therefore confirms once more the above general analysis results. 
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Appendix 



The formalism of the fundamental multivaluedness, providing true quantum 
chaos in Hamiltonian systems, is based on the division of an "unresolvable" 
problem into two parts. One of them, represented by the modified Schrodinger 
equation, eq. (5), is relatively simply structured, easily tractable and gives 
finally the fundamental dynamic uncertainty in the form of problem splitting 
into plural realisations. The second part is represented by the auxiliary system 
of equations, eqs. (4), which is more difficult for the detailed analysis (for 
example, because of the large number of the participating equations), but 
influences the first part and the solution of a problem in a less direct, largely 
qualitative way. It provides a necessary 'material support' for the fundamental 
uncertainty and plays a minor (though irreducible) role for the basic origins of 
dynamical randomness. Nonetheless, one needs to know the general properties of 
the auxiliary system solutions as well as the effective approximations that 
could be used in practical applications of the formalism of FMDF. We 
demonstrate these features of the auxiliary system of equations for the more 
transparent time-dependent case, while the time-independent problem can be 
analysed in the same way and with the same results accompanied by minor 
technical changes. 

A useful approximation for the auxiliary system, eq. (4b), can be obtained 
if we rewrite it in the form: 

Mra)V^J(ra) + n(ra)V^J(ra) = eWkiro) - V^r(ra)V^L'(ra) (A.1) 

k^^ k 

(we recall that by definition k, k' ^ 0). We can apply now the same EP formalism 
that was used in section 2.1 to obtain the modified Schrodinger equation, 
regarding the sum in the right-hand side of eq. (A.1) as a perturbation. The result 
differs from that of eq. (5) in two points: the EP in this case evidently depends 
on k, and it cannot be obtained in a closed form, but only as a series.*) Omitting 
the detailed calculations quite similar to those of section 2.1, we write down 
the first term of this series for the EP kernel: 

Vr(ra)V-r(ra1V^o„(ra)<*(ra1 

. (A.2) 

^k-k' - ^on + h(On(k - k') 

n 



Vyfc(ra,ra') = 




*) The latter feature seems to have a fundamental implication involved with the fractal structure of 
dynamic complexity. The dependence of the higher-order terms on the eigenvalue to be found accounts 
for the appearance of the respective higher-order fractal branches, together with similar dependence of 
another quantity, the 'averaged energy addition' (see section 4. Ill for more details). However, these 
details should not directly influence the structure of the current level of the structure of chaos. 
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where {el„} and {y/ln) a''© the eigenvalues and the eigenfunctions for the 
unperturbed potential: 

ho{ra)y/lnira) + Vo{ra)vLira) = £^¥^(^0) ■ (A.3) 
As £^ = £a - hcOnk, one can rewrite eq. (A. 2) in the form: 




Vk{ra,ra) = Z y ] ■ (A.4) 

n 

This leads to the following graphical-analysis representation (see eqs. (13), 
(14)) of the corresponding modified Schrodinger equations: 

,,,. 0,0 00 /A i;\ 



where 




k'i^k 

^nn{k,£kn) — ^^^^^ o , , ' (A-6) 



-'k n 



and 



= J fi?raV^o«'(»*a)^'fc(»*cj)V^fc„(ra) , 

So 

£kn = |^i?raV^^*(ra)[/lo(ra) + Vo(ra)] V/'L(ra) 



The singular denominators in eq. (A.6) determine the asymptote positions 
(section 2.2), and we see that, contrary to the case of eqs. (13), (14), now there 
is no additional splitting of solutions: one has A^a asymptotes for each k which 
provide only one realisation (at th\s level of complexity, cf. section 4. III). 
Moreover, the obtained eigenvalue set depends little on k (it is rather shifted as 
a whole with changing k). The weak dependence on k of the corresponding EP can 
be also traced from eq. (A.4). If one adds another terms of the series for the EP, 
this will not change our qualitative results because each of these terms 
contains the same denominator as eq. (A.6), but raised to the corresponding 
power, and the respective additional summations of terms quadratic in Vk^^cs) in 
the numerator. Note, by the way, that this situation gives an idea about how one 
can reproduce regularity, within our method, for an arbitrary integrable system. 
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Another approximation and a way of analysis of tlie auxiliary system can 
be obtained if we add the same quantity, y/l{rc) XI ^k'i^a), to both sides of the 

initial eq. (A.1): 

ho{ra)¥kira) + \Vo{ra) + V^J(ra) = 

(A.7) 

= e>;j(ra) - llvkira)[Wk-k^{ra) - Wki^o)] ■ 

k^^ k 

If we suppose that the effective potential, 

Vi(ra) - Vo(ra) + L ^kiro) = V{r^,t=0) - Vk{r^) , (A.8) 

k^^k 

as well as the solution, \j/l{ra), depend only weakly on k, this will provide a self- 
consistent approximation for eq. (A.7): 

ho{ra)wl{rc) + Vi(ra)V^J(ra) = e>J(ra) . (A.9) 

This approximation is evidently more precise for the 6-like perturbation 
dependence on t containing many sufficiently large higher harmonics (in the 
limit it coincides with the 5-functional kicks of the standard model). Contrary 
to this, the previous series expansion of the EP is more relevant to the case of 
quasi-harmonic perturbation. However, the qualitative properties of the 
auxiliary system solutions are similar for both approximations. The difference 
is determined by the effective potential which changes approximately between 
Vo(ra) for a harmonic perturbation and Vi(ra) for 6-like kicks. These two 
limiting EP models show weak dependence on k and differ mostly by their 
amplitudes: AVi » AVq. In this sense one may say that the EP follows, to some 
degree, the form of the perturbation. 

In the time-independent case similar analysis leads to analogous results 
concerning the qualitative properties of the auxiliary system solutions. We do 
not reproduce it here noting only that the detailed expressions deal with 
rather than So (see eqs. (3)). 
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